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We study a type la supernova explosion using three-dimensional numerical simulations based on 
reactive fluid dynamics. We consider a delayed-detonation model that assumes a deflagration-to- 
detonation transition. In contrast to the pure deflagration model, the delayed-detonation model 
releases enough energy to account for a healthy explosion, and does not leave carbon, oxygen, and 
intermediate-mass elements in central parts of a white dwarf. This removes the key disagreement 
between simulations and observations, and makes a delayed detonation the mostly likely mechanism 
for type la supernovae. 



PACS numbers: 97.60.Bw, 26.30.+k, 47.70.Fw, 47.40.-x 

Type la supernovae (SN la) are produced by ther- 
monuclear explosions of white dwarf (WD) stars com- 
posed primarily of C and O nuclei and detached degener- 
ate electrons. The most probable explosion scenario in- 
volves a binary star system, in which a WD can increase 
its own mass by accreting material from its compan- 
ion until it approaches the Chandrasekhar limit, 1.4M0. 
Near this limit, a small increase in mass results in sub- 
stantial contraction and compression of the WD. The 
compression increases the temperature, accelerates nu- 
clear fusion reactions, and triggers the thermonuclear 
runaway [1] that eventually ignites thermonuclear burn- 
ing near the WD center. This starts a thermonuclear 
explosion that releases ~10^^ ergs during a few seconds. 
The energy is produced by a network of thermonuclear 
reactions that begins with the original ^^C and ^^O nuclei 
and ends in the formation of ^^Ni and other iron-group 
elements. Considerable amounts of intermediate-mass el- 
ements (IME), such as Ne, Mg, Si, S, and Ca, are created 
as well. Type la and other types of supernovae play an 
important role in stellar nucleosynthesis and in releasing 
the newly synthesized elements into interstellar medium, 
thus providing raw material for next generations of stars 
and planets. 

Because of their extreme and predictable luminosity, 
SN la are extensively used as standard candles to measure 
distances and estimate cosmological parameters critical 
for our understanding of the global evolution of the uni- 
verse. To improve these measurements, we need compre- 
hensive theoretical and numerical models of SN la that 
describe details of the explosion and connect them to ob- 
served characteristics of SN la, such as spectra and light 
curves. One-dimensional (ID) numerical models have 
been extensively used to test general ideas about possi- 
ble explosion mechanisms [2-7]. In particular, ID models 
have ruled out the possibility of a thermonuclear det- 
onation, a supersonic shock-induced combustion mode, 
as a sole mechanism for SN la explosions. A detonation 
propagating through a high-density WD produces mostly 
Ni and almost none of IME observed in SN la spectra. 
One-dimensional models have also shown that a detona- 



tion can produce IME if it propagates through a low- 
density WD preexpanded during the initial deflagration 
stage of the explosion. These delayed-detonation models 
[8-14], which have a deflagration-to-detonation transition 
(DDT) at some stage of the thermonuclear explosion, are 
the most successful in reproducing observed characteris- 
tics of SNe la. Many crucial physical details, however, 
including the mechanism of DDT and the turbulent flame 
structure, are missing by definition from ID models be- 
cause SN la explosions are intrinsically three-dimensional 
(3D) phenomena. 

Full-scale 3D numerical simulations of thermonuclear 
supernova explosions have become a reality during the 
past few years [15-18]. They have shown, in particular, 
that the development of a turbulent thermonuclear flame 
in the gravitational field of a WD allows funnels of un- 
burned and partially burned material to remain in the 
vicinity of the WD center until the end of the explosion. 
This would produce distinct signatures of low-velocity 
C, O, and IME in SN la spectra. As the observed spec- 
tra do not show these signatures, the deflagration model 
must be incomplete. Previously we suggested [15,18] but 
did not prove that a detonation triggered by the turbu- 
lent flame could burn the remaining material near the 
WD center and make the model consistent with observa- 
tions. Here, we test this hypothesis using a 3D numerical 
delayed-detonation model of SN la explosion in which a 
deflagration undergoes a transition to a detonation. 

Input Physics and Numerical Implementation. — The 
numerical model discused in details in [15,18] is based on 
reactive Euler equations that include gravity terms and 
are coupled with an equation of state for a degenerate 
matter and a simplified kinetics of energy release. The 
equations are integrated on a Cartesian adaptive mesh 
using an explicit, second-order, Godunov-type numeri- 
cal scheme. The model describes compressible fluid dy- 
namics on large scales in an exploding WD including the 
expansion of a star, Rayleigh- Taylor (RT) and Kelvin- 
Helmholtz (KH) instabilities, turbulence, pressure waves, 
shocks, and detonations. The nuclear kinetics is approxi- 
mated by a four-equation mechanism [8,15] that describes 



the energy release, consumption of C, nuclear statistical 
quasi-equilibrium (NSQE or QSE) and nuclear statisti- 
cal equilibrium (NSE) relaxations, and neutronization. 
The turbulent flame speed is provided by an additional 
subgrid model [15,18] that takes into account physical 
processes at scales smaller than the computational cell 
size. In particular, it assumes that turbulent burning on 
small unresolved scales is driven by the gravity-induced 
RT instability. 

The model is able to reproduce the two different 
regimes of the thermonuclear burning in a WD, a sub- 
sonic deflagration and a supersonic detonation. These 
regimes differ by the mechanism of propagation of the 
reaction wave: a deflagration involves heat conduction or 
turbulent mixing, and a detonation involves shock com- 
pression. For both regimes, the energy is released by the 
same network of thermonuclear reactions, and the phys- 
ical thickness of the reaction front strongly depends on 
density. It can be up to 12 orders of magnitude less than 
the WD radius Rwd for deflagrations [19,20] and up to 
10 orders of magnitude less than Rwd for detonations 
[21,22]. Since the large-scale simulations described here 
do not resolve length scales smaller than \Q^^Rwd, the 
reaction fronts at high densities are still unresolved. We 
explicitly resolve only parts of the reaction zone associ- 
ated with NSE relaxation that become very large at low 
densities and cause an incomplete burning that produces 
Si and other IME. Resolution tests show that the mini- 
mum computational cell size dx^m = 2.6 x 10^ cm used 
here for the deflagration stage and dxmin = 10.4 x 10^ cm 
used for the detonation stage are adequate for this type 
of simulations. 

Deflagration Stage. — The initial conditions for the 
deflagration stage model a IAMq WD in hydrostatic 
equilibrium with initial radius Rwd = 2 x 10^ cm, ini- 
tial central density pc = 2 x 10^ g/cm^, spatially uni- 
form initial temperature T = 10^ K, and uniform initial 
composition with equal mass fractions of ^^C and ^^O 
nuclei. The burning was initiated at the center of WD 
by filling a small spherical region at r < Q.QlbRwD with 
hot reaction products without disturbing the hydrostatic 
equilibrium. We model one octant of the WD assuming 
mirror symmetry along the cc = 0, y = and z = 
planes. The computational domain is a cube with a side 



of Xr. 



5.35 X 10» cm. 



The development of the thermonuclear flame was de- 
scribed in detail earlier [18]. The initially spherical flame 
propagates outwards with a laminar velocity ^^ 100 kni/s, 
becomes distorted due to the RT instability, and forms 
multiple plumes at different scales. Buoyancy causes the 
hot, burned, low-density material inside the flame plumes 
to rise towards the WD surface. The same gravitational 
forces also pull the cold, unburned, high-density material 
between the plumes down towards the center. The flame 
becomes turbulent and forms a dynamic convoluted sur- 
face penetrated in all directions by very irregular funnels 
of unburned and partially burned material. 

The higly developed 3D flame surface increases the 



burning rate and improves the estimation of energy re- 
lease compared to ID and 2D deflagration models. The 
intense convection on large scales also causes the burned 
material to spend less time in high-density central parts 
of WD, thus reducing the neutronization [15]. The same 
convective flows, however, bring unburned material to the 
central parts of the exploding star. As a result, substan- 
tial amounts of C, O, and IME remain near the WD cen- 
ter by the end of the explosion. As we have shown [15,18], 
this makes predictions of the 3D deflagration model in- 
consistent with observed spectra of SN la. 

Detonation Stage. — The disagreement between pre- 
dictions from the pure deflagration simulation and ob- 
servations strongly suggests that the turbulent flame in 
SN la triggers a detonation. The process of DDT involves 
events occurring at small scales that are comparable to 
the detonation wave thickness, and, thus, cannot be di- 
rectly modeled in large-scale simulations. To study the 
effects of a detonation, we therefore assume a time and 
a location for DDT. (A similar approach has been used 
previously in ID [8-13] and 2D [23-25] delayed-detonation 
models.) Now the deflagration results are initial condi- 
tions, and we impose a hot spot to ignite the detonation. 
The time and location for the detonation initiation are 
parameters that can be varied and optimized. Here, we 
explore the three cases (a), (b), and (c) defined below. 

The case (a) corresponds to central detonation initi- 
ation at 1.62 s after the beginning of the deflagration. 
By that time, 1/3 of WD mass has burned, the WD 
radius has increased by a factor of 1.55, and the den- 
sity of unburnt material near the center has dropped to 
2.5 X 10^ g/cm^. A detonation at this density produces 
mostly Ni and propagates outwards at '--^12, 000 km/s, 
which is comparable to the expansion velocities induced 
by subsonic burning. When the detonation reaches un- 
burned material with density below (1 — 5) x lO'' g/cm'^, it 
begins to produce IME. Different parts of the detonation 
front that exit different funnels collide with each other, 
coalesce, and eventually reach the surface of the star. 

The detonation transforms all C and O in central parts 
of the WD into iron-group elements, and produces IME 
in outer layers. This drastically changes the distribution 
of nuclei compared to that produced by the pure defla- 
gration. Funnels of unburned C and O disappear from 
central parts of the WD. Iron-group elements form a dis- 
tinct core surrounded by a layer of IME. Angle-averaged 
mass fractions of the main elements calculated for the 
deflagration and the delayed-detonation models are com- 
pared in Fig. n 

Similar results were obtained for the delayed- 
detonation case (b), for which the detonation was initi- 
ated at 1.62 s at 10* cm off center and produced a moder- 
ate asymmetry in composition. The results indicate that, 
during the period of detonation propagation, the density 
of the expanding unreacted material ahead of the shock 
can decrease by an order of magnitude compared to its 
value at the end of the deflagration stage. Because the 
detonation burns material to different products at dif- 
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FIG. 1: Angle-averaged mass fractions of the main elements 
as functions of scaled distance from the WD center produced 
by the deflagration (d) and delayed-detonation (a) models 
at 1.94 s after the beginning of the explosion. The delayed- 
detonation model corresponds to case (a) described in the 
text. Xmax ~ 5.35 X 10* cm. Lines marked as Mg, Si, and 
Ni represent estimated cumulative mass fractions of elements 
from Ne to Mg, Si to Ca, and Ti to Ni, respectively. The esti- 
mations are based on a four-equation nuclear kinetic scheme 
[15,18] and the reaction zone structure of a ID detonation 
wave calculated in [21] with a detailed nuclear kinetics. 



ferent densities, it can create a large-scale asymmetry in 
composition if it starts far from the WD center. A similar 
conclusion based on 2D simulations was made in [25]. 

The asymmetry effect in our simulations is limited be- 
cause we calculate only one octant of a WD and impose 
mirror boundary conditions. The degree of asymmetry 
would increase if the simulations were performed for a full 
star. Then the second mirror-reflected spot for detona- 
tion initiation would be eliminated. Three-dimensional 
simulations [15,26] also show that a developing flame, 
unrestricted by mirror boundaries, can move away from 
the WD center, thus creating a large-scale asymmetry at 
very early stages of the explosion. 

For the case (c), the detonation was initiated at the 
WD center at 1.51 s when 1/4 of WD mass has burned, 
the WD radius has increased by a factor of 1.30, and the 
density of unburnt material near the center has decreased 
to 4.4 X 10® g/cm'^. Case (c) produced more iron-group 
elements than cases (a) and (b) because the detonation 
propagated through higher-density material. The earlier 
detonation initiation also resulted in a faster explosion 
that released 15% more energy. Total energies for all 
three cases and the deflagration model are shown in Fig. El 
as functions of time. The total energy Etot here is the 




FIG. 2: Total energy as function of time for deflagration (d) 
and delayed-detonation cases (a), (b), and (c) described in 
the text. Energy units are 10^" ergs. 



difference between the energy released by thermonuclear 
reactions and the binding energy of the star. Eventually 
Etot will be transformed into kinetic energy of expanding 
material that can be measured in observations of SN la. 

Figure El shows that the total energy released by 
delayed-detonation models, (1.3 — 1.6) x 10^^ ergs, is 
much higher than the energy released by the deflagration 
model ~0.6 x 10^^ ergs. The reason for this is that the 
deflagration is able to burn only about a half of the WD 
mass. The rest of the material expands to the densities 
below ~10^ g/cm^ that do not support the thermonu- 
clear burning. A detonation propagates faster and burns 
almost all of the WD mass before the material expands to 
low densities. The total energy released by the delayed- 
detonation models is in agreement with a typical range 
(1 — 1.5) X 10^^ ergs obtained from SN la observations 
[7]. 

Discussion and Conclusions. — Figure Q shows the 
averaged distribution of elements at 1.94 s, the time 
when the WD surface reached the computational domain 
boundary, but the detonation did not yet reach the WD 
surface. At this time, the detonation front propagates 
through low-density outer layers of the star and produces 
mostly IME. All iron-group elements have already formed 
at higher densities. The total mass of iron-group ele- 
ments created by the explosion is 0.78, 0.73, and 0.94 
solar masses (Mq) for delayed-detonation cases (a), (b), 
and (c), respectively. Most of this mass is the radioac- 
tive ^^Ni that provides the energy source for the observed 
luminosity of SN la. The mass of '^^Ni estimated from 
observational data is about O.6Af0 for a typical SN la 
[27], and is in agreement with the total mass of iron- 
group elements produced by delayed-detonation models. 
For the deflagration model, the total mass of iron-group 
elements is only 0.47 Mq, which is insufficient to account 
for the luminosity of a typical SN la. 

The carbon-oxygen layer that remains between the 
detonation front and the WD surface will continue to 
burn as the detonation advances. Oxygen in outer lay- 
ers, which expand to densities below ~10^ g/cm'^ be- 



fore the detonation reaches them, will remain unburned. 
Carbon is likely to remain unburned for densities below 
(1 — 3) X 10^ g/cm'^. Unburned C and O in outer lay- 
ers would produce spectral signatures only in the high- 
velocity range. 

There are possibilities, however, for a delayed detona- 
tion to leave small amounts of C and O in inner parts of 
WD. For example, a detonation propagating through a 
thin, sinuous funnel of unburned material can fail if the 
funnel makes a sharp turn. A developing turbulent flame 
can also disconnect some funnels from the rest of the 
unburned material, thus creating unburned pockets that 
cannot be directly reached by a detonation wave. These 
pockets may or may not ignite when strong shocks gener- 
ated by detonations reach them. The cellular structure of 
thermonuclear detonations in carbon-oxygen matter [22], 
and the ability of cellular detonations to form pockets of 
unburned material that extend far behind the ID reac- 
tion zone, can also contribute into incomplete burning. 
All these phenomena occur at length scales comparable 
to the reaction zone thickness that are not resolved in 
large-scale simulations reported here, and require addi- 
tional studies. 

There have recently been efforts to detect low- velocity 
C in SN la spectra that could result from the funnels 
of unburned material near the WD center [28] . The re- 
sults [28] indicate that C can be present at velocities 
11 000 km/s. Even though this velocity is much lower 
than 20 000-30 000 km/s usually attributed to C in SN la 
spectra [29-31], it is still too high for the ejecta formed 
from central parts of a WD. For C and O, spectral sig- 
natures are difficult to observe, and estimated veloci- 
ties of these elements are subject to large uncertainties. 
Intermediate-mass elements, however, produce distinct 
spectral lines and their velocities are well defined. The 
minimum observed velocities for IME [28,32] are large 
enough {r-^ 10 000 km/s for Si) to rule out the presence 
of IME near the WD center, as is predicted by the defla- 
gration model. A discussion on this subject can also be 
found in the recent article [33]. 

Figure n shows that, in contrast to the 3D deflagra- 
tion model, the 3D delayed-detonation model of SN la 
explosion does not leave C, O, and IME in central parts 
of a WD. This removes the key disagreement between 
simulations and observations, and makes the 3D delayed 
detonation a promising mechanism for SN la explosion. 
Further analysis of 3D delayed detonations on large scale 
requires 3D radiation transport simulations to produce 
spectra, and a detailed comparison between the calcu- 
lated and observed spectra of SN la for different initia- 
tion times and locations. The uncertainty in detonation 
initiation can only be eliminated by solving the DDT 
problem that involves physical processes at small scales. 
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